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Abstract 

We discuss synchronization in networks of neuronal oscillators which are in- 
terconnected via diffusive coupling, i.e. linearly coupled via gap junctions. In 
particular, we present sufficient conditions for synchronization in these net- 
works using the theory of semi-passive and passive systems. We show that the 
conductance-based neuronal models of Hodgkin-Huxley, Morris-Lecar, and the 
popular reduced models of FitzHugh-Nagumo and Hindmarsh-Rose all satisfy 
a semi-passivity property, i.e. that is the state trajectories of such a model 
remain oscillatory but bounded provided that the supplied (electrical) energy is 
bounded. As a result, for a wide range of coupling configurations, networks of 
these oscillators are guaranteed to possess ultimately bounded solutions. More- 
over, we demonstrate that when the coupling is strong enough the oscillators 
become synchronized. Our theoretical conclusions are confirmed by computer 
simulations with coupled Hindmarsh-Rose and Morris-Lecar oscillators. Finally 
we discuss possible "instabilities" in networks of oscillators induced by the dif- 
fusive coupling. 
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1. Introduction 



Synchronous behavior is witnessed in a variety of biological systems. Exam- 
ples include the simultaneous flashing of fireflies and crickets that are chirping 
in unison |39j , the synchronous activity of pacemaker cells in the heart [2a ] and 
synchronized bursts of individual pancreatic /3-cells 



2G 



331 ] . For more examples see 



381 ] and the references therein. It is well known that individual neurons 



in parts of the brain discharge their action potentials in synchrony. In fact, 
synchronous oscillations of neurons have been reported in the olfactory bulb, 
the visual cortex, the hippocampus and in the motor cortex g, |34| . Presence 



or absence of synchrony in the brain is often linked to specific brain function or 
critical physiological state (e.g. epilepsy). Hence, understanding conditions that 
will lead to such behavior, exploring the possibilities to manipulate these condi- 
tions, and describe them rigorously is vital for further progress in neuroscience 
and related branches of physics. 

We present results on synchronization of ensembles of neuronal oscillators 
which are being interconnected via gap-junctions, i.e. a linear electrical cou- 
pling of the form g ■ (Vi(t) — V^i)) where the constant g represents the synaptic 
conductance and V\ (t) — V<x (f ) denotes the difference in membrane potential of 
the neurons at the pre-synaptic side and the post-synaptic side at time t, respec- 
tively. Recently it has been pointed out that gap-junctions play an important 
role in synchronization of individual neurons [2]. 

Several attempts have been made to understand when synchronization of 
neurons coupled via gap junctions occurs. In P, 14, 16, [2 42 1 phase equa- 



tions and phase response curves are used to analyse neurons coupled via gap 
junctions. They all conclude that for increasing coupling the synchronous state 
becomes stable. However, the use of phase equations is only justified when the 
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coupling between the cells is weak. In general, the results for strong coupling 
are rare J]. In J] Coombes uses a piecewise linear model of spiking neurons 
which allows to extend the results for weak coupling (using phase equations) for 
strong coupling. Chow and Kopell [3j used Integrate-and-Fire kind of models 
to investigate synchronization via gap junctions. They showed using spike re- 
sponse functions (for the Integrate-and-Fire models an analytic expression for 
this function exists) that, depending on the shape of the spikes, the firing fre- 
quency and the coupling strength, stable phase locked states exist. When the 
coupling is large the oscillators will synchronize. They showed using simulations 
that for more realistic models the results hold true as well, however no rigorous 
mathematical proof is presented. In [131 ] conditions for synchrony in two coupled 
Hodgkin- Huxley neurons are presented; If the coupling between the neurons is 
strong enough, then the neurons will synchronize. In 22 1 synchronization for 



multiple interconnected chaotic Hindmarsh-Rose neurons is discussed. 

We will generalize the results obtained in [22j and present conditions for syn- 
chrony of diffusively coupled identical neuronal oscillators for general network 
topology. From the zoo of models of neuronal activity (see [1 11 ] for a review) we 
will focuss on four popular oscillators, namel y th e conductance based, biophysi- 
cally meaningful models of Hodgkin-Huxley [kJ and Morris-Lecar [18| , and the 
more abstract models derived by FitzHugh-Nagumo and Hindmarsh-Rose 

[9]. First we demonstrate that, despite the difference in the range of behavior 
that these models are capable to produce, these models have an important col- 
lective property. This property is that each model is semi-passivqj. Second, 
using the concept of semi-passivity, introduced in 30(, we will show that a set 



of these diffusively coupled neuronal oscillators will always possess bounded so- 
lutions. Next, under condition that the coupling between the neurons is large 



'we will formally introduce semi- passivity in Definition 12. II in Section [2] 
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enough, i.e. there is a high-conductive pathway between the neurons, we show 
that the oscillators will become synchronized. 

This paper is organized as follows. In Section [3 we introduce the notion 
of semi-passivity and we show that four models mentioned above are all semi- 
passive. Next, in Section [31 a theorem adopted from [28[ is presented which 



provides sufficient conditions under which the oscillators show synchronous be- 
havior. We demonstrate in Section [3.21 using computer simulations that ensem- 
bles of Hindmarsh-Rose and Morris-Lecar oscillators will end up in synchrony 
whenever the coupling between the neurons is large enough. In Section [3] we 
briefly discuss that it is not obvious that systems being interconnected via dif- 
fusive coupling will have bounded solutions and eventually end up in synchrony. 
In particular, we show that two "dead" cells can become "alive" when being 
interconnected via diffusive coupling, i.e. the cells start to oscillate due to the 
interaction. Finally, Section [5] concludes of the paper. 

Throughout this paper we use the following notations. The symbol R stands, 
as usual, for the real numbers, R+ denotes the following subset of R: R+ = 
{x G R|x > 0}. The Euclidian norm in R" is denoted by |j-||, ||x|| 2 = x T x 
where the symbol T stands for transposition. The symbol I n defines the n x n 
identity matrix and the notation col (x\, . . . , x n ) stands for the column vector 
containing the elements xi, . . . , x n . A function V : R™ — + R+ is called positive 
definite if V(x) > for all x & R n \ {0}. It is radially unbounded if V(x) — > oo 
if ||x|| — > oo. If the quadratic form x T Px with a symmetric matrix P = P T is 
positive definite, then the matrix P is positive definite, denoted as P > 0. The 
symbol C r denotes the space of functions that are at least r times diffcrcntiablc. 
Consider k interconnected systems and let Xj denote the state of a single system, 
then the systems are called synchronized if limj_ >00 — = 0, i,j € 

{1,2,...,*}. 
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2. Semi-passivity 

We represent a neuronal oscillator as the general system 



x = f(x) + Bu, 

(1) 

y = Cx, 



where state x £ R n , input m £ ft is an depolarizing or hyperpolarizing (input) 
current and output y G R denotes the membrane potential of the neuron. Fur- 
thermore, / : R™ — > R™ is a C 1 -smooth vector field and the matrices B and C 
are of appropriate dimensions. 



Definition 2.1 (Passivity and semi-passivity [43, 28]). The system fl} is 
called 

i) passive in T> C R" if there exists a nonnegative function V : T> — > R + , T> 
is open, connected and invariant under the dynamics (jX|) 7 V(0) — 0, such 
that the following dissipation inequality 

V(x) = ^^(f(x)+Bu)<y T u (2) 
holds; if V — R™ the system is called passive; 



ii) semi-passive in T> if there exists a nonnegative function V : T> C R™ — > 
R+, T> is open, connected and invariant under V(0) = 0, such that 

V(x) = ^^(f(x)+Bu)<y T u^H(x), (3) 

where the function H : T> C R™ — > R is nonnegative outside the ball B 
with radius p 

3p>0, \\x\\>p=>H(x) > g(\\x\\), 

with some nonnegative continuous function g(-) defined for all \\x\\ > p; if 
T> = R™ the system is called semi-passive; 



Hi) strictly semi-passive (in T>) if the function H(-) is positive outside some 
ball BcV. 

A semi-passive system behaves similar to a passive system for large enough 
|x||. Hence a semi-passive system that is interconnected by a feedback u = <p(y) 



satisfying y T <p(y) < has ultimately bounded solutions [43 



i.e. regardless 



how the initial conditions are chosen, every solution of the closed-loop system 
enters a compact set in a finite time and stays there, see Figure [TJ Moreover, 
this compact set does not depend on the choice of initial conditions. 
Consider k identical neuronal oscillators of the form 



%i = f( x j) + Bu ^ 



(4) 



where j = 1, . . . , k denotes the number of each system in the network, Xj E R" 
the state, Uj E R the input and yj E R the output of the j th system, i.e. 
the membrane potential, smooth vector field / : R n — > R n and vectors B — 
[10 ... 0] T and C — [1 ... 0] are of appropriate dimensions. Note that many 
neuronal models are in this form or can be put in this form via a well-defined 
change of coordinates. 

The k neurons ^ are coupled via diffusive coupling, i.e. a mutual intercon- 
nection through linear output coupling of the form 



Ijk (Vj - Vk) 



(5) 



where jji — jij > represents the synaptic conductance and — yj is the 
difference in membrane potential of neurons i and j. 
Defining the k x k coupling matrix as 



r = 



Eh 
3=2 Tlj 



"712 



"721 



-7fei 



Ek 
7 = 1,7^2 72j 



"7fc2 



-Ilk 
-Ilk 



Efe-1 
, = 1 Ikj 



(G) 
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the diffusive coupling functions (|5|) can be written as 



u = -Ty (7) 

where u = col (ui, . . . , Uk) and y — col (yi, . . . , j/fe). Since T — T T all its eigen- 
values are real and T is singular because all rowsums equal zero. Moreover, 
applying Gerschgorin's theorem (cf. js?]]) about the localization of the eigen- 
values, it is easy to verify that T is positive semi-definite. We assume that the 
network cannot be divided into two or more disconnected networks. Hence the 
matrix T has a simple zero eigenvalue. 

Proposition 2.1. Consider a network of k diffusively coupled systems (|1|) 7 (|5|). 

Assume that each system in the network is semi-passive, then the solutions of 
all connected systems in the network are ultimately bounded. 



Proof The proof is adopted from |28j. Let the j th system in the network be 
semi-passive with the storage function V(xj), where Xj is the state of the j th 
system. Denote W{x) = J2j=i V( x j) where x — col (xi, . . . , Xk), then 

W(x) = ]T V( Xj ) < ]T yJ Uj - H( Xj ) - -y T Ty - ]T H( Xj ) < 0, (8) 

outside some ball in R nfc . Note that the quadratic term y T Ty is nonnegative 
since T is semi-positive definite. This directly implies that the solutions of the 
interconnected systems are bounded and exist for all t > to. 

Remark 2.1. Even if the systems are not identical, but each individual system 
is semi-passive, then the network will still have bounded solutions. This follows 
directly from i.e. the storage function for the network is simply the sum of 
the storage functions of the individual oscillators. 

Remark 2.2. Consider a collection of k neurons that interact via chemical 
synapses, where the chemical synapse is modeled as the Fast Threshold Modula- 
tion (FTM) coupling introduced in \3aj . i.e. 

k 

Ui = -9ijH{yj -0)(a- y l ), (9) 
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where gij £ R>o denotes the synaptic conductance, aGtis the synaptic rever- 
sal potential which determines whether the synapse is inhibitory or excitatory, 
and the function H(-) is typically chosen as the Heaviside function such that 
neuron j will influence neuron i only if the membrane potential of neuron j ex- 
ceeds some threshold 9 6 R. It is not hard to verify that semi-passive neuronal 
oscillators interconnected via chemical synapses have bounded solutions. (This 
follows from the fact that "Y^yiUi < outside some ball in R fc ; i.e. the "supplied 
energy" is bounded.) 

We are now ready to prove that the neuronal models of Hodgkin-Huxley, 
Morris-Lecar, FitzHugh-Nagumo and Hindmarsh-Rose all satisfy the semi-passive 
property. Hence the solutions of networks of these oscillators with a diffusive 
coupling exist and are bounded. 

Hodgkin-Huxley model 

The most important model in computational neuroscience is probably the 
Hodgkin-Huxley model [lOj . Consider the Hodgkin-Huxley equations : 

Cxi = g Na xlx 3 (E Na - xt) + g K x\ (E K - xi) + g L (E L - x\) + I + u 



±i = ai{x\) (1 - Xi) - f3i(xi)xi, i — 2,3,4 

with y = x\ is the membrane potential, state iC^C R 4 , input u G R, positive 
constants gNa, 9k, 9l, C € R and constants /, E^ a , Ek, El 6 R. The functions 
ctj(-) and are defined as 



(10) 



25 - s 




oti{s) 



«3(s) 



0.07e- s / 20 , 
10 - s 



100 ( e ( 1 - s / 1 °) - l) 



(11) 



ft 00 



ft 00 



1 



ft(s) 



e (3-s/10) + I 

0.125e" s/8 °. 
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The states Xi represent so-called activation particles which satisfy Xi(t) £ (0, 1) 
for all t > to whenever Xi(to) £ (0, 1). 

Proposition 2.2. The Hodgkin- Huxley model is semi-passive in X where 

X = {x £ R 4 |0 < Xi < 1, i = 2, 3, 4}. (12) 

Proof First, we will prove that for all to < t\, to, t\ £ R: 

CI) X\{t) exists on the interval t £ [to, t\] and remains bounded if the input u 
is bounded; 

C2) Xi{t) £ (0, 1) on the interval t £ [t , h] provided Xi(t ) £ (0, 1). 

We do so by invoking a contradiction argument. Suppose that CI) does not 
hold. Let us denote 

u* = sup |K*)||. (13) 

te[to,*i] 

According to assumptions of the proposition such u* must exist. The right-hand 
side of (fl0|) is locally Lipschitz, hence its solutions are defined over a finite time 
interval. Let [to, T] be the maximal interval of their existence. Let us pick some 
arbitrarily large constant M £ R+. Then there should exist a time instant t[ 
such that 

\\x{t)\\>M, Vt>t[. (14) 
Consider the internal dynamics 

±i = a>i(xi) (1 - — Pi{xi)xi, i = 2,3,4. (15) 

One can easily verify that ai(x\) > 0, (3i{x\) > for all (bounded) x\. Hence 
on the boundary Xi = we have ii > and at the boundary xi = 1 we have 
Xi < 0, i.e. Xi can not cross the boundaries. Hence the set (0, 1) is forward 



9 



invariant under the Xi dynamics, i.e. for all Xi(to) G (0, 1), 



< asi(*) < 1, Vte[i ,T]. (16) 
Then, according to (HHJ), (flt?|) the following holds 



<e- A ( t - t «)|x 1 (t )|+p+^*, Vte[t ,T] (17) 



where p, A are positive constants of which the value do not depend on M. 
Combining (TTi|) and (fT7)l we obtain 



M< \\x(t)\\ <e- x ^\ Xl (t )\+p+ju*, Vieft.T] (18) 



where M is arbitrarily large and p, xi(to), and 1/Xu* are fixed and bounded. 
Hence we have reached contradiction, and CI) hold. This automatically implies 
that C2) holds too. 

To finalize the proof of semi-passivity of (|10p , consider the storage function 
V : X -> R+, V = ^x\ + | E xl Then 

, , 4 , (19) 

+ {gNaX 2 xzE Na + g K x±E K + g L E L + I) x x 

4 

- ^ fai^i) - \) 2 - jj + (3i{xi)xi 2 
Note that — fa,(xi) ^(a;, — i) 2 — |J + /^(iei)^ 2 ^ < for each x t outside (0, 1). 
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Because C2) holds we obtain 



V <x\u — ghx\ + c\X\ 



4 



(20) 



i=2 



where constant 



Cx = max \dig N aE N a + d 2 g K E K + 9lE l + 1\ X 



di,d 2 G[0,l] 



(21) 



x sign (dig Na E Na + d 2 g K E K + ghE L + I) . 



Given that (|2"0)) holds for all t, the Hodgkin-Huxley model is semi-passive in X. 

Morris-Lecar model 

The Morris-Lecar model is a planar system that models the voltage 
oscillations in the barnacle giant muscle fiber. The Morris-Lecar model is given 
by the following equations 

C±x = g L (E L -x%) + gcaaoo (xi) (E Ca -x±)+ g K x 2 (E K -xx)+I + u, 
X 2 =1] (xi) (Poo(xi) - x 2 ) , 



with y = x\ denoting the membrane potential, state x £ X C K 2 , input «el, 
constant parameters El, Ec a , Ek £ H, positive constants gL, gca, gK £ K. 
and functions 



(22) 




(23) 
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with fj > 0, Ei, E%, E3, E4, fj £ R. Like in the Hodgkin-Huxley equations, the 
states xi represent an activation particle which satisfies x 2 {t) £ (0, 1) for all 
t > to provided X2(to) £ (0, 1). 

Proposition 2.3. The Morris-Lecar model is semi-passive in X where 

X = {x £ R 2 |0 < x 2 < 1}. (24) 

Proof Notice the forward invariance of the set (0, 1) under the ^-dynamics. 
The proof is similar to the proof for the Hodgkin-Huxley equations. 



FitzHugh-Nagumo model 

The FitzHugh-Nagumo model 



19j is one of the simplest models of the 



spiking dynamics of a neuron. The model is given by the following set of differ- 
ential equations 

x 3 

Xl = Xi — - — x 2 + 1 + u, 

3 (25) 

±2 = 4> {xi + a - bx 2 ) , 

where y = Xi represents the membrane potential, state x = (xi, x 2 ) T £ R 2 , 
input u £ R and positive constants a,b, <p £ R. Constant parameter I 6 1 
determines the output-mode of the model (either spiking or quiet). 

Proposition 2.4. The FitzHugh-Nagumo equations satisfy the semi-passivity 
property ([3]). 

Proof Consider the storage function V : R 2 — > R+ 

V=\{4 + \4). (26) 

Then 

T 4 

V = xiu — + Xi + Ixi — bx 2 + ax 2 - (27) 
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Therefore V (x 1,0:2) < x\u— H(x\, X2) with H(xi, x%) = ^■—x\—Ixi J \-bx\~ax2, 
i.e. the FitzHugh-Nagumo neuron is semi-passive. 



Hindmarsh-Rose model 



Consider the Hindmarsh-Rosc 



.1 



equations 



x\ = —ax\ + bx\ + X2 — X3, + I + u 

±2 = c — dx\ - x 2 (28) 
±3 = r (s (xi +w) — x 3 ) 

where y — x\ represents the membrane potential, state x = (x\, X2, x$) T £ R 3 , 
input «£l and constant positive parameters a, 6, c, d, r, s, w 6 R. The constant 
parameter / 6R determines again the output-mode of the model, which in this 
case, depending on the choice of parameters, can be resting, bursting or spiking. 
Moreover, for some parameters it can even behave chaotically. 
Proposition 2.5. The Hindmarsh-Rose model is semi-passive. 



Proof The proof is adopted from 22 1 . Consider the storage function V : R 3 — > 
R+ 

V=\{xl + »xl + ±xl) (29) 
with constant /i > 0. Hence 

V = x\u — ax\ + bx\ + X1X2 + Ix\ + ficx2 — [idx\x2 — \ix\ + wx% — ~x$. (30) 

Let 

- ax\ - \1dx\x2 = -aXtxf - o(l - Ai) (x\ + 2a( ^ Xl) x 2 ^j + 4^1- (31) 
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and 



- ilX% + XXXI = -\l\ix\ - - A 2 ) (X2 - 2ai (i_a 2 ) x i ) + 4»(l-\ 2 ) x l ( 32 ) 

with Aj G (0, l)cE,i = l, 2. Then 
V —x±u 

- a\ x x\ + bx\ + 4fi ^_x 2 ) xl + Ixi 

- (/ iA 2 - 4a(l-A 1 ) ) X 2 + MCX 2 

V 7 (33) 

- - s x\ + wx 3 

- ^(1 - A 2 ) (x 2 - a^ib^ i) 
-a(l-A 1 )(^ + 5ZT ^ I77a;2 ) 2 . 

Let /i < 4aA 2(4- A i) . Then it follows directly that the Hindmarsh-Rose model 
satisfies the semi-passivity property ([3]). 

Remark 2.3. Many biophysically meaningful neuronal models, i.e. conduc- 
tance based models like the Hodqkin- Huxley and Morris-Lecar models, share the 
same structure, see for instance fjj}, [Tj, IJuj - In particular, the evolution of the 
membrane potential is given by an equation of the form 

k 

Cv(t) = u(t)+J2lj(t) (34) 

3=1 

where n£E denotes the membrane potential, C € R>o is the membrane capac- 
ity, u G R is the input and ionic currents Ij(t) = gj(t)(Ej — v(t)) with constant 
reversal potential Ej G R and time-varying conductance Qj(t) > for all t. The 
conductance is typically given as 

m 

'/, </..lR (35) 

i=l 

with maximal conductance cjj G R>o, nonnegative integers pij and voltage de- 
pendent gating variables Si(v(t)), where the gating variables satisfy Sj(t) G (0, 1) 
for all t > t whenever Si(t ) G (0, 1). 

All models of neuronal oscillators of this form are semi-passive fin E x 
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(0, 1) x ... x (0, 1)), and the proof for semi-passivity is similar to the proof 
presented for the Hodgkin- Huxley model. 

Remark 2.4. Consider the class of Integrate- and- Fire neurons, i.e. neuronal 
models of the form 

x = a — bx + u, if x > xthres, then x <— c, (36) 

where the output y — x represents the membrane potential, u is the input, pos- 
itive constants a, b, Xthres is the threshold potential and c is the value to which 
the membrane potential x is reset to after firing. The state of an Integrate- 
and-Fire neuron will always be bounded, i.e. c < x < x t hres, hence we do not 
need a semi-passivity argument to guarantee the solutions of such a model to be 
bounded. 



3. Synchronization of diffusively coupled neuronal oscillators 

In the previous Section we showed that the solutions of diffusive coupled 
neurons (of the Hodgkin-Huxley, the Morris-Lecar, the FitzHugh-Nagumo and 
the Hindmarsh-Rose type) remain bounded. Using these results we provide 
conditions for which the neurons end up in synchrony. 

Since the matrix CB is nonsingular, the systems ([4]) can be transformed into 
the following form 

Vj = a(Vj>%j) + CBuj = a(yj,Zj) + u j} 

(37) 

where yj G R, Uj G R, Zj G R m , m = n — 1, and sufficiently smooth functions 
a : R x R m -> R, q : R m x R -> R"\ 

Theorem 1. \2a] Consider the k systems (|37p and assume that: 
i. each system 

Vi = a(yj,Zj) + Uj, 
z, } = q{z h yj), 

is strictly semi-passive; 
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ii. there exists a C 2 -smooth positive definite function Vq : R m — » R+ and a 
positive number a £ R suc/i t/ia£ f/ie following inequality is satisfied 

(VV Q (z' - z")) T (q(z>, y>) - q(z", y')) < -a \\z> - z"\f (39) 

for all z', z" G R m cmd j/' £ R. 

Then, for all positive semi-definite matrices V all solutions of the closed-loop 
system (|37|) . ([7| are ultimately bounded. Let the eigenvalues Xj of T be ordered 
as = Ai < A2 < • • ■ < Afe. TTien ttere exists a positive number A suc/i t/iat i/ 
A2 > A t/iere exists a globally asymptotically stable subset of the diagonal set 

A = {ijj £ R, Zj e R m : Vi = yj,Zi = zj, i, j = 1, . . . , k} . (40) 

Remark 3.1. One can easily verify that Theorem]]] remains true in case that 
each system {57J is semi-passive in T>, for T> as defined in Definition \2.1\ 

According to Theorem [T] the problem of examining the asymptotic stability 
of the synchronized state of all oscillators in the network is reduced to 

i. verification of the assumptions for an individual oscillator, and 

ii. computation of the eigenvalues of the coupling matrix T. 

It follows that if for a given network topology the coupling is large enough, i.e. 
A2 exceeds the threshold A, then the neurons will synchronize. Moreover, once 
the threshold value A is known one can easily determine whether the neurons in 
networks with different topologies synchronize or not by computing the eigen- 
values of the corresponding coupling matrix. This is the Wu-Chua conjecture 



44|. The effect of the network topology on the synchronization can also be 
investigated using, for instance, the Connecting Graph Stability method [lj]. 

3.1. Convergent systems 

There exists a sufficient condition to check whether inequality (|39[) of The- 
orem [T] is satisfied or not. Therefore, let us introduce the notion of convergent 
systems. 

Definition 3.1 (Convergent systems). @, [H/ Consider the system 

z = q{z,w(t)), (41) 
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where the external signal w(t) is taking values from a compact set W C R. The 
system (|41|) is called convergent if 

i. all solutions z(t) are well-defined for all t G (—00, +00) and all initial 
conditions z(0), 

ii. there exists an unique globally asymptotically stable solution z w (t) on the 
interval t € (—00, +00) from which it follows 

lim \\z(t) - z w (t)\\ =0 (42) 

t — >oo 

for all initial conditions. 

The long term motion of such systems is solely determined by the driving input 
w(t) and not by initial conditions z(0), i.e. the systems "forget" their initial 
conditions. A sufficient condition for a system to be convergent is presented in 
the next lemma. 

Lemma 1. U If there exists a positive definite symmetric m x m matrix 
P such that all eigenvalues Xi(Q) of the symmetric matrix 

Q(z,w) = ^ 

are negative and separated from zero, i.e. there is a S > such that 

Xi(Q)<-5<0, (44) 
with i = 1, . . . , m for all z € R m , w £ W, then the system (|41[) is convergent. 

It follows that if there exists such a matrix P such that each system ij = q{zj , yj) 
satisfies (|4"5]) . (|31)l. i.e. each system Zj = q(zj,yj) is convergent, then inequality 
(f39|) of Theorem [T] is satisfied. 

One can easily verify that the internal dynamics of the models of Hodgkin- 
Huxley, Morris-Lecar, FitzHugh-Nagumo and Hindmarsh-Rose are convergent, 
(use P = I in (|43f and the result follows.) 

3.2. Illustrative examples 

In the previous section we have shown that all four the models satisfy the 
semi-passivity condition. Moreover, the internal dynamics of these systems are 



P 



dq 



[z,w] 



Oq 



(z,w] 



P 



(43) 
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equivalent to a convergent system. Therefore, according to Theorem[TJa network 
consisting of the presented oscillators shows bounded solutions and, in case 
the coupling is strong enough, all oscillators will end up in perfect synchrony. 
However, the goal is here not to determine the exact threshold values for which 
the network starts to synchronize. Such threshold values can be expressed in 
terms of the system parameters (see, for instance 



ICS 

.11 



22j or [l| for Hindmarsh-Rose 



neurons), or they can be determined by computing, for instance, the transversal 
Lyapunov exponents of the coupled systems [24J. Here, the goal is only to show 
that for large enough coupling the neurons will synchronize. 

Synchronization of Hindmarsh-Rose oscillators 

Consider a network of eight diffusively coupled Hindmarsh-Rose neurons 

= —aXj tl + bx 2 x + Xjfi - Xj^ + I + uj 
xj,2 = c — dx 2 3 l - Xj :2 (45) 
Xj,3 =r(s (xj.i +w) - x j<3 ) 

where j — 1, ... ,8 denotes the number of the oscillator in the network. We use 
the following set of parameters: a = 1, 6 = 3, c = 1, d = 5, r = 0.005, s = 
4, w = 1.6180, I = 3.25. With these parameters each Hindmarsh-Rose neuron 
has chaotic solutions [9(. Let the eight oscillators be connected as shown in 



Figure 2(a) with corresponding coupling matrix 
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(46) 



The smallest nonzero eigenvalue of T 1 is A3, « 2.587. O ur simulations show 
that the neurons synchronize when 7 > 0.387, which corresponds to A 1 = 1.00. 
(This agrees with the numerical results obtained in, for instance, Q] , where it is 
shown that two diffusively coupled Hindmarsh-Rose neurons synchronize when 
the coupling strength 7 > 0.50, i.e. A = 1.00.) Figure [3] shows the simulation 
results of the network of Hindmarsh-Rose oscillators with coupling 7 = 0.39 such 
that A2 ~ 1.01. The top panel shows the x\ states of the eight oscillators, the 
middle panel shows the X2 states and the X3 states are depicted in the bottom 
panel. The first 500 [s] the systems are uncoupled and one sees the systems are 
not synchronized. After 500 [s] the coupling becomes active, indicated by the 
arrows in Figure [31 and all systems rapidly synchronize. 

Synchronization of Morris-Lecar oscillators 

Next we synchronize eight Morris-Lecar oscillators which are connected ac- 
cording to the graph depicted in Figure[2(b)| The corresponding coupling matrix 
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is given as 



p 2 = 
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such that the smallest nonzero eigenvalue of T 2 is X\ w I.277. Each Morris- 
Lecar oscillator is given by the following set of equations 



Cx jt i =g L (E L - xjs) + gcaaoc (xj,i) (E Ca ~ Xj,i) + 

+ QKXj.i{E K - xjs) + I + Uj, (48) 

with j = 1, . . . , k denoting the number of the oscillator in the network and 
functions 



1 ( , fs-Ex 

= 2i 1+tanh h^ 



aoo(s) 

0oo(s) 

rj(s) = fj cosh 



1 ( , fs-E 3 

= 2 1+tai HV 



E. 



(49) 



2E 4 



We used C = 1, ,g L = 0.5, £ L = -50, g Ca = 1-1, £ C « = 100, g K = 2, 
= -50, I = 30, Ex = -1, E 2 = 15, £3 = 0, £4 = 30, 77 = 5 in our numerical 
simulations. Figure [4] shows the simulation results for the eight diffusively cou- 
pled Morris-Lecar oscillators with 7 = 0.01. The first 250 [s] the oscillators are 
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uncoupled and do not synchronize. Then, after 250 [s] the coupling is turned on, 
which is again indicated by the arrow, and all oscillators become synchronized. 



4. Diffusion driven instabilities 

In this section we show using two simple examples that it is not trivial that 
systems interacting via diffusive coupling have bounded solutions and possibly 
end up in synchrony. In particular, we demonstrate that diffusive coupling 1) 
can make the solutions of the interconnected systems to become unbounded, 
and 2) can make systems, which have an asymptotically stable equilibrium in 
isolation, to produce stable oscillations. 

Example 4.1 (Unbounded solutions). Consider the linear (non-minimum 
phasqjj) stable transfer function 



H(s) 



1 



2s 2 + 2s + 1 



A possible state space realization for the system is 



where 



A = 



X 


= Ax 


+- B 


u, 


y = 


Cx, 
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Consider now two diffusively coupled systems ([51]) 

i\ — Axi + -fBC(x 2 - xi), 
x 2 = Ax 2 + jBC(xi - x 2 ), 



(50) 

(51) 
(52) 

(53) 



Clearly the origin of each uncoupled system is globally asymptotically stable. 
However, the system is not semi-passive, and when 7 > 0.65 12 [ for 7 = 0.6512 
the system undergoes a Poincare-Andronov-Hopf bifurcation the solutions 
of the interconnected systems become unbounded. 



2 a system is non-minimum phase if it has unstable zero dynamics, i.e. the internal dynamics 
with constraint y(t) = are unstable, cf [27j |. 
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Example 14.11 shows how the diffusive coupling between two not semi-passive 
systems results in unbounded solutions. A similar phenomena is encountered 
in networks of diffusively coupled Chua circuits, cf. [4l| . The piecewise linear 
model of the Chua circuit is not semi-passive (the Chua attractor is not globally 
stable) and due to the interaction the trajectories of the systems can be driven 
outside the domain of attraction such that the solutions grow unbounded. 



The following example is taken from [27j. It shows how two systems, which 

both have an asymptotically stable equilibrium in absence of interaction, start 

to produce stable oscillations when the systems interact via diffusive coupling. 

Example 4.2 (Diffusion driven oscillations). Consider two systems which 
interact via diffusive coupling: 



where matrices A, B and C are as presented above. Again the origin of an 
isolated system is asymptotically stable and when 7 = 0.6512 the (linearized) 
system undergoes a Poincare-Andronov-Hopf bifurcation. Hence the coupled sys- 
tems (|54p start to produce stable oscillations whenever 7 > 0.6512, see Figured 
for simulation results. 

The key mechanism for the oscillations is the Poincare-Andronov-Hopf bifurca- 
tion and the non-minimum phaseness of the systems. The diffusive interaction 
between initially silent cells is essential for generating stable oscillatory behav- 
ior in some neuronal (and other biological) systems, see [15I and the references 
therein. The authors demonstrate that the main reason for the oscillations is 
that the internal variables, e.g. (in) activation particles, have the tendency to 
oscillate. However, these oscillations are being suppressed through a negative 
feedback mechanism. The diffusive coupling will destroy the feedback mecha- 
nism causes the internal variables and, hence, the membrane potential to start 
to oscillate. Note that the mechanism is the same as in our example, i.e. the 
internal dynamics are not minimum phase. However, the goal here is not to 




±1 = Ati(1 + ||iz;i|| 2 )+7BC(x2-:ei), 
x 2 = Ax 2 (l + \\x 2 \\ 2 ) + jBC(x! - x 2 ), 



(54) 
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discuss the machinery for the generation of these oscillations in detail. We refer 
the reader to 



27| for more details. 
Note that the four models described above do have minimum phase internal 
dynamics since the internal dynamics are convergent. Hence no "spontaneous" 
oscillations due to diffusive interaction will occur in networks of Hodgkin-Huxley, 
Morris-Lecar, FitzHugh-Nagumo and Hindmarsh-Rose neurons. 



5. Discussion 

We have presented sufficient conditions for synchronization in networks of 
diffusively coupled neuronal oscillators. The results are constructive in the fol- 
lowing sense: 

i. we have considered different classes of neuronal oscillators, i.e. neuronal 
oscillators whose dynamics are described via the Hodgkin-Huxley formal- 
ism and neuronal oscillators of the FitzHugh-Nagumo and Hindmarsh-Rose 
type; 

ii. we have shown that all these oscillators are semi-passive with a quadratic 
storage function. A consequence is that, when semi-passive systems are 
being coupled via coupling of the type ([5]), the network possesses bounded 
solutions (see Proposition ^. ip . 

iii. the internal dynamics of all these neuronal oscillators are convergent, i.e. 
the subsystem z = q(z, y) of each oscillator satisfies the conditions as stated 
in Lemma [T] As a result (|3"5|) of Theorem [TJ is satisfied; 

iv. since the oscillators are semi-passive and the internal dynamics are con- 
vergent it is possible, according to Theorem [TJ that all oscillators in the 
network end up in stable synchrony. The criteria for synchronization is 
that for a given network topology the strength of the interconnections is 
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large enough since topology and coupling strength influence the smallest 
nonzero eigenvalue of the coupling matrix T. 

Theorem Q] allows to decompose the problem of finding (sufficient) conditions 
for stable synchronization in the network of k coupled oscillators into some con- 
ditions of the individual oscillators (semi-passivity and internal convergent dy- 
namics) and conditions on the network (topology and coupling strength, which 
both influence the smallest nonzero eigenvalue of the coupling matrix). Our 
results can therefore be applied to neuronal networks interconnected via strong 
coupling and with general network topology. Our theory supports the result of 
[13I and the simulation results of Chow and Kopell 3] for strong coupling (note 
that the model of the interneuron and the Traub- Miles model are both semi- 
passive and have convergent internal dynamics). Moreover, our results hold even 
when the neurons behave chaotically. In [6| the authors discuss the emergence 
of clusters in all-to-all coupled networks as function of the coupling strength. 
Those clusters might emerge when the coupling is not strong enough to end up 
in synchrony. The emergence of clusters for diffusively coupled neurons satisfy- 
ing the assumptions of Theoremf! can be explained for general network topology 



using the theory discussed in 



29. 



3 11 ] - The goal of this paper is to show that 



neurons interconnected via gap junctions will posses bounded solutions and, 
moreover, synchronize whenever the coupling strength is large enough. Deter- 
mining sharp synchronization thresholds however will still depend highly on the 
type of neurons involved and their specific set of parameters. We have also 
shown that diffusively coupled systems which are not semi-passive might have 
unbounded solutions. A probably more interesting property, at least from the 
biological point of view, is that diffusively coupled non-minimum phase systems 
which are initially silent can start to produce stable oscillations. 
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In this paper we considered networks of diffusively coupled neurons without 
any time-delay. However, in a physical system one would expect that it takes 
some (small amount of) time to transmit a signal. Therefore it is interesting 
to analyse synchronization in networks where time-delays are included in the 
coupling. Sufficient conditions for synchronization in time-delayed networks are 
presented in, for instance, [21|. Here, the semi-passivity property in combina- 
tion with a small-gain theorem provides a sufficient condition for boundedness of 
the trajectories of the coupled systems. Next sufficient conditions for synchro- 
nization in terms of Linear Matrix Inequalities (LMIs) are derived. However, 
solving the LMIs is computationally involving, especially when the networks be- 
come large and complicated. Moreover, the results might be very conservative. 
It would be interesting as well to investigate the emergence of stable synchro- 
nization in pulse-coupled networks. This is because most neurons are actually 
coupled via so-called chemical synapses, i.e. an impulsive type of coupling. It 
is shown in \v\ that certain Integrate-and-Fire neurons in a network with all- 
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3G 



32] 



to-all connections synchronize for almost any initial condition. In 
synchronization of more realistic neuronal oscillators in pulse-coupled networks 
is discussed. However, rigorous constructive results about what conditions the 
oscillators should satisfy and the effect of a particular network topology on the 
synchronization are not present nowadays. Even for systems that can be rep - 



resented as the seemingly simple (pulse-) coupled Kuramoto oscillators, cf 
the problem of global synchronization is not tackled in full generality. Networks 
with strong interactions and/or chaotic regimes remain problematic. It is for 
future research to explore the possibilities in these topics. 
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Figure 1: Semi-passivity; every solution enters the ball ||x|| < p in finite time and stays there 
as time increases. 




(a) Graph 1 




(b) Graph 2 

Figure 2: Eight diffusively coupled oscillators. Each interconnection has weight 7. 
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xi = fi(x 1 ,x c ) + u 



X2 = fa{X2,x{) 



Xn fn [Xn ; X \ ) 



